Effect of Hund coupling in the one-dimensional SU(4) Hubbard model 
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The one-dimensional SU(4) Hubbard model perturbed by Hund coupling is studied, away from 
half-filling, by means of renormalization group and bosonization methods. A spectral gap is always 
present in the spin-orbital sector irrespective of the magnitude of the Coulomb repulsion. We 
further distinguish between two qualitatively different regimes. At small Hund coupling, we find 
that the symmetry of the system is dynamically enlarged to SU(4) at low energy with the result 
of coherent spin-orbital excitations. When the charge sector is not gapped, a superconducting 
instability is shown to exist. At large Hund coupling, the symmetry is no longer enlarged to SU(4) 
and the excitations in the spin sector become incoherent. Furthermore, the superconductivity can 
be suppressed in favor of the conventional charge density wave state. 
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I. INTRODUCTION 

The interplay of spin and orbital degrees of free- 
dom plays an important role in diverse correlated elec- 
tron systems.^ Recently, the one-dimensional (ID) spin- 
orbital models have been studied intensively motivated 
by the discovery of the quasi- ID spin-gapped materi- 
als, Na2Ti2Sb20 and Na2V205i^ These materials can be 
modeled by a quarter-filled two-hand Hubbard model. 
Even in this approximation, the situation is rather com- 
plex owing to the large number of independent coupling 
constants in the problem. 

A first attempt to understand the effect of the high de- 
generacy induced by the orbital degrees of freedom was 
to consider the situation where both orbital and spin de- 
grees of freedom play an identical role. In this case, the 
two-band Hubbard model possesses a large symmetry and 
it becomes SU(4) symmetric in the spin-orbital sector 
and depends on only one coupling: the Coulomb repul- 
sion U . At quarter-i\\\ing it was found that spin-orbital 
and charge degrees of freedom separate at low energy. 
The spin-orbital sector remains massless for all values of 
U > Q and displays quasi long-range antiferromagnetic 
order with three critical modes whose dynamics is de- 
scribed by a SU(4) level 1 (fc=l) Wess-Zumino-Witten 
(WZW) model. This result is in agreement with the ex- 
act result by Sutherland for the SU(4) Heisenberg chain-. 
What happens for the U(l) charge excitations strongly 
depends on U . It was found that a Mott-Hubbard tran- 
sition from a massless phase at small U to an insulating 
phase at large U takes place at a nonvanishing critical 
value Uc = 2.8t4i Of course, the SU(4) symmetry in the 
spin-orbital sector is not likely to be present in real mate- 
rials, so that a systematic study of the effects of possible 
symmetry breaking operators is necessary to account for 



the experimental results eventually. 

One of the simplest symmetry breaking is to break 
the SU(4) symmetry down to SU(2)spin x SU(2)orbitai- A 
detailed renormalization group studj*^ revealed that the 
SU(4) symmetry is dynamically enlarged at low energy. 
Furthermore, the massless phase in the spin-orbital sector 
survives in an extended region of coupling constant space. 
These results were confirmed numericallyji^ 

In the present work we investigate a different symme- 
try breaking perturbation which is always present in real 
materials, the Hund coupling. We also consider the cases 
of general fillings away from the half-filled one. In com- 
parison with the spin-orbital model studied in Refs.0|^, 
the Hund coupling breaks the SU(4) symmetry further 
down to SU(2)spinXU(l)orbitai- As for our most impor- 
tant result, we find that a spectral gap opens in the spin- 
orbital sector for an arbitrarily small Hund perturbation. 
We further distinguish between two qualitatively differ- 
ent regimes. At small Hund coupling, we find that the 
SU(4) symmetry is dynamically enlarged at low energy 
like in the spin-orbital model. We further show that a 
superconducting (SC) instability is present in the charge 
sector. At large enough Hund perturbation, the SU(4) 
symmetry is no longer fully enlarged. Instead, we find 
a partially SU(2) symmetry enlargement in the orbital 
sector. In this phase, the SC instability may disappear 
in favor of the conventional charge density wave (CDW) 
instability. 

This paper is organized as follows. In Sec. II, we 
present our model and discuss its symmetry properties. 
The renormalization group analysis at weak coupling is 
performed in Sec. HI, where we also discuss the physi- 
cal properties in both spin-orbital and charge sector. In 
Sec. IV, we develop a complementary strong coupling 
approach for the quarter-Q\\eA case. We conclude this 
paper in Sec. V. 
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II. MODEL HAMILTONIAN AND ITS 
SYMMETRY PROPERTIES 

The Hamiltonian we consider is the U(4) Hubbard 
model with a Hund couphngr^ 



H = Ha + Hj, 



(1) 



with 



Hn = 



E (-C+icLc.+io'.+H.c.) 



u 



and 



-2JE S» 



(3) 



where the Ciaa are the electron operators at the site i 
in orbital a — (1,2) with spin a. In Eq. SiQ, = 

Scrcr' cLa [f Ciaa', dcnotc thc Spin 1/2 operators of 
electrons in both bands a — (1,2). We further assume 
that the hopping is diagonal in orbital space, i.e. if f+i = 
tSaa', and that U and J are positive. 

The symmetry properties of Eq. are most clearly 
seen by introducing the U(l) charge Q and the SU(4) 
spin-orbital generators T"^, A= {!,■ ■ ■ , 15), as follows: 

Q ~ E] '^lacr'^iacri 

= E i.'.'[MXf c.,a^. (4) 

i,a' a' acr 

where are the generators of SU(4) Lie algebra^"^. A 
convenient explicit realization of thc is 

i=(^)^'52', ^^:'(y)a', n/2(^)-'(^)S', (5) 

where cr° and t"", a = 1,2, 3, are the Pauli matrices acting 
on spin and orbital space, respectively. An appropriate 
labeling of the SU(4) generators in Eq. ((SJ is as follows. 
To each SU(4) index ^ = (1, • • • ,15) we associate a pair 
of indices, such that (a, b) ^ (0, 0), (a, 6 = 0, 1, 2, 3) with 
the convention that cr^ = r*^ = Id2. For an example, 
the first three generators of Eq. ((SJ can be alternatively 
expressed as M^'^'^\ 

When J = 0, using Eqs. (|4I5|I . the Hamiltonian H is 
clearly seen to commute with both Q and all of T"*'s, 
thus it is U(4)=U(l)chargcXSU(4)spin_oibitai Symmetric. 

The Hund coupling Hj does not affect the charge sec- 
tor but breaks the SU(4)spin-orbitai symmetry. Indeed, 
when J 7^ 0, in addition to the obvious SU(2) invari- 
ance in spin space generated by M(°'°\ (a — 1,2,3), 
H is also invariant under the U(l)orbitai group in or- 
bital space generated by M^*^'^\ Thus, the Hund cou- 
pling breaks the SU(4)spin-oibitai symmetry down to 



SU(2)spinXU(l)orbitai- In Comparison with the spin- 
orbital model studied in Refs. 5 6] and Ref.F^I, the Hund 
term breaks the symmetry SU(2)orbitai further down to 

U(l)orbital- 



III. RENORMALIZATION GROUP ANALYSIS 
AT WEAK COUPLING 

The effective low energy theory associated with the 
Hamiltonian ^ is obtained by taking the continuum 
limit in a standard way. At small enough (U/t, J/t) and 
at low energy, the electron operators Ciaa may be ex- 
panded near the two Fermi points ±kp: 



(6) 



where oq is the lattice constant. In the continuum limit, 
the effective Hamiltonian is expressed naturally in terms 
of the chiral U(l) and SU(4) current densities4i 



J, 



L{R) 



E^I 

E ^lim,.'.'[M-^ro:fMn),c..- (7) 



a' a' a<j 



After some algebra, discarding irrelevant operators, we 
find that the charge degrees of freedom decouples from 
the spin-orbital ones away from half-filling: 



where 



and 



He = dx 



H = He + Hs 



TTV. 



-fiJR + Ji)+9cJRJL 



(8) 



(9) 



Hso — 



AiE[-^r^r 



aO) 
a3) 



JaQ) jiaO) 
■^R '^R 



■^R '^R 



a 

5^E^ 

a 

hY.W'j'R^''+jt^''j'R^ 

a 

hY.W'j'R 

a 

~ \ .(0,1) .(0,1 

54 \Jl Jr 



33) 



1) , .(0,2) .(0,2) 

•^L 'Jr. 



~ .(0,3) .(0,3) 
^ 95 J L "Jr 



(10) 
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In Eqs. H9I10|I . Vc = wf(1 + 3Uao/2TrvF) and Vso = 
vf{1 — Uao/2TrvF) are the charge and spin-orbital veloc- 
ities, where vp = 2tao siiikpao is the Fermi velocity. We 
observe that all interactions of Eqs. (|9I1U|) are marginal 
and of the current-current type.^*' Therefore, the low en- 
ergy physics will result from a delicate balance among 
the different interaction terms in Eq. (|10|l . The bare 
coupling constants in both charge and spin-orbital sec- 
tors are expressed in term of U and J as follows: 

3 

g^ = -Uao, (11) 



and 



Ai = -A2 = Jao, gi = 2UaQ + 2Jaa, 
g2 = 2Uao + Jao, gs = 2Uaa - 2Jao, 
54 = 2?7ao - 3Jao, gz^2Uan. (12) 

The effective Hamiltonian in the charge sector [Eq. 0] 
is that of Luttinger liquid: 



2 



dx 



K, 



(13) 



where 4>c — 4>cL + 0c_r and 9^ = 4>cL — (t>cR are the charge 
boson field and its dual field, respectively. The chiral 
boson fields 4>c,l/r are defined in terms of currents as 
follows: 



(14) 



Jl{R){x) = V - d^(t)cL{R){x). 



Therefore, the charge sector is massless and the low 
energy properties are determined by the nonuniversal 
charge exponent Kc which is given, at leading order in 
f/, by: 



3Uao 



TIVp 



-1/2 



< 1. 



(15) 



The charge velocity can be reexpressed as — vp/K"^. 
The situation at hand is similar to what happens in the 
SU(4) Hubbard model at quarter filling. 

The effective Hamiltonian in the spin-orbital sec- 
tor is that of the SU(4)i Wess-Zumino-Novikov-Witten 
(WZWN) model with the central charge c = 3, perturbed 
by marginal interactions. This is similar to what happens 
in the spin-orbital model studied in Refs.'s"^ and Rcf."?!. 
Due to the complexity of the interaction pattern, namely 
the five coupling constants instead of three in the spin- 
orbital model, the situation in the spin-orbital sector in 
the presence of a Hund term requires a careful analysis of 
the renormalization group (RG) flow. Out of the seven 
coupHng constants entering in Eq. ifTUI) . the Ai and A2 
terms are purely chiral and are not renormalized at lead- 
ing order. Furthermore, they do not influence the scaling 
of the cjiS. One- loop RG equations are easily found by 



current algebra method>ii*i2i 

dgi 

dt 
dg2 

dt 
dgs 

dt 
dgi 

dt 
dg5 

dt 



2 n 2 2 

-91 - 2^2 - 93, 



-2gi.g2 - 9295 - 5354, 
-2gig3 - 25234, 
-3g2.g3 - 9i95, 



-^91 



where 



9i 



9I, 



9i 



(16) 



(17) 



and t = \nL \s the RG time. 

We have performed a detailed numerical analysis of the 
RG flow associated with Eqs. \i<6\ . In the following, we 
summarize our results. 

When J = 0, the interaction is irrelevant for [/ > 
and all coupling constants flow toward the SU(4)i fixed 
point at g* = 0. There are no other fixed points associ- 
ated with Eq. H16|l . One of our most important results 
is that a nonzero value of the Hund coupling, J 7^ 0, 
destabihzes the SU(4)i fixed point and drives the system 
toward strong coupling. This indicates that a gap opens 
in both spin and orbital sector with Mgpin ~ -Morbitai ~ 
exp(— C/J), where C is a positive constant of order t. 
The present situation is completely different from the 
one encountered in the spin-orbital model where the crit- 
ical SU(4)i phase was not entirely destabilized by the 
SU(2)xSU(2) symmetry breaking perturbation.^ Though 
a gap opens in the spectrum, the low energy effective the- 
ory still depends on the relative magnitude, i] — J /2U , 
between the Coulomb repulsion and the Hund coupling. 
Indeed, one finds two qualitatively different behaviors of 
the RG flow depending on 77. 

-The regime A: The SU(4) symmetric regime. The 
first is a regime with the enlarged SU(4) symmetry. This 
regime occurs for 77 ^ 1 or [/ ^ J. This regime will be 
referred to as the regime A from now on. 

With these initial conditions, though all the coupling 
constants gi{t) — > ±cxd when t t* , they asymptotically 
match the following particular RG-invariant ray: 



91 = -92 = -93 = 54 = 9a 



(18) 



On that ray, omitting chiral terms, one may write Eq. 
(|10|l in an explicit SU(4) invariant form: 



dx 



[2ttv*. 



~* jA jA 

9 •Jl "Jr 



(19) 

where we have performed the duality transformation for 
a, 5^ 0: 



a.h 



a.h 



ja,b 



(20) 
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and V* is an effective spin-orbital velocity. In fact, there 
is a velocity anisotropy in the model. We find, how- 
ever, that to leading order in J/t such an anisotropy of 
velocities scales to zero at low enough energyji^ Thus, 
we find that the symmetry is dynamically enlarged to 
SU(4) to the one-loop accuracy. Of course the validity 
this result which relies on the loop expansion may be 
questioned j^ii^ii^ but it is reasonable to conjecture that 
the enlargement of the symmetry is likely to hold beyond 
the perturbation theory. In any case such an (enlarged) 
SU(4) symmetry is meant to be approximate in the sense 
that small corrections to pure SU(4) behavior should be 
expected due to the neglected irrelevant operators. 

-The regime B: The SU(2)orbitai enlarged regime. The 
second is a regime B where SU(2)orbitai symmetry is par- 
tially enlarged [from U(l)]. This regime occurs for the 
large Hund coupling J 3> U or rj ^ 1. For 77 ^ 1 we 
find that the SU(4) symmetry is no longer fully enlarged, 
and instead we observe a partial SU(2)orbitai symmetry 
enlargement in the orbital sector. With the initial condi- 
tions satisfying U <^ J, the RG flows drive the coupling 
constants to a regime where: 



-00 



gi{t) < g^it) « .92(i) < gsit) « gi{t) < 0. (21) 

In this regime the RG equations Eq. H16(l can be ap- 
proximately decoupled. Indeed, at long RG time, the 
contributions of 32 and 93 to the RG equation for (74 and 
55 can be neglected and one obtains: 



dg4_ 
dt 

dg5 
dt 



-9495, 



2 

^9a, 



(22) 



which are nothing but the RG equations of the U(l) 
Thirring model in the orbital sector with effective initial 
conditions |54(i)| > |.95(^)|- In this regime, it is known 
that the SU(2) symmetry is restored at larger RG time. 
Once the anisotropy between 34 and g^ becomes small, 
so does the anisotropy between 52 and g^ as can be seen 
from the equation 



d{g2 - .93) 
dt 



94(92 - 53), 



(23) 



since (32 — 53) < and 34 < 0. Therefore, in the strong 
coupling regime, the effective Hamiltonian approximately 
depends on three independent coupling constants: 



[54 = 55] < [92 = .93] < gi < 0. 



(24) 



With the above relation Eq. (|24|l . the interacting part 
of Eq. (irU)) displays an SU(2)spinX SU(2)orbitai symme- 
try. It is important to notice that no further symmetry 
restoration is expected since the coupling constant in the 
orbital sector 54 is much larger than the one in the spin 
sector gi . This behavior is in contrast with what happens 
in the A phase for 77 ^ 1. 

The study of the physical properties near the bound- 
ary between the A and B regimes is a nontrivial problem. 



Whether they are separated by a quantum phase transi- 
tion point or they are smoothly connected by a crossover 
region can answered only by methods far beyond the per- 
turbation theory. This problem will be addressed else- 
where. Within the one-loop accuracy, we find that the 
RG flow qualitatively changes from the A-type to the 
B-type as 77 decreases below rjQ ~ 0.5, which is a reason- 
able value, but we were not able to conclude in favor of 
a quantum phase transition. 



A. Physical Properties and Order Parameters 

1. Spin- Orbital Sector 

In order to get some insights in the physical properties 
of the low energy physics it is appropriate to change the 
parametrizations of the fluctuating fields. We first notice 
that the two sets of SU(4)i currents, J^'°^ and jjj^'^j, a = 
(1, 2, 3), span two spin and orbital SU(2) algebras. More 
precisely, they are SU(2)fc=2 currents. This stems from 
the fact that the SU(4)fc=i WZW model is equivalent to 
the sum of two decoupled SU(2)fc=2 WZW modeliSiii As 
in Ref.0 we shall take advantage of the representation 
of the SU {2)k=2 algebra in term of three (real) Majorana 
fermionsiSiiS: 



J 



(a,0) 



2 = 



J' 



(a,b) 
R{L) 



^ ^abc cb 

2 

* abc (b 

2 ?t,fl.(L) 'St,-R(i)' 



^CsML)i'tML)AaM^Q. (25) 



where ^^"^^^^^ and £,1 



t,R{L)^ 



- (1, 2, 3), are the Majorana 
fermions associated with the spin and orbital degrees of 
freedom. In term of these Majorana fermions the effective 
theories in both A and B regimes take a nice form. 

In the regime A, the effective low energy Hamiltonian 
can be obtained from Eq. H1U|) with the condition 91 — 



-92 



-.93 = .94 = .95 = .9 < imposed: 



n 



a 

a 

5[$^«-<) 



(26) 



where k^^^^ 



is{t)Ris{t)L- "^^^ Hamiltonian is easily 
seen to be S0(6) invariant upon a duality transformation 
in the orbital sector: 



it 



t,R 



(27) 



which is the equivalent of Eq. (|20|) when Eq. H2U|) is 
expressed in terms of the Majorana fermions. Under 
the duality transformation Eq. (|27|l , the Hamiltonian 
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becomes that of the integrable SO (6) Gross-Neveu 
(GN) modeliSiSi which has been first obtained in Ref. @ 
as the effective low energy theory for the massive phase 
of the spin-orbital model. In this respect we find that, 
though the Hund coupling Eq. (PJ breaks the original 
SU(4) symmetry further than spin-orbital model, it is 
not responsible for the new low energy physics as far as 
J is not too large in the spin-orbital sector. Therefore, 
many of the conclusions drawn in Ref.|^|^|3| still hold 
for moderate values of J. In particular, the excitation 
spectrum consists of a kink and an anti-kink with mass 
m along with a fundamental fermion of mass V2m. The 
existence of the fundamental fermion as a stable quasi- 
particle implies that the spin excitations are coherent: 
a sharp resonance at w = v^ti is expected to appear 
in the dynamical structure factor, in particular, the 2kp 
component of spin-spin correlation function. This can be 
checked by explicit calculation via order/disorder opera- 
tor formalism of Ising model [Eqs. (|43l I44|l ] . 

In the regime B, where J ^ U, the effective Hamilto- 
nian is that of two coupled SO (3) GN models, one in the 
spin sector and the other in orbital sector: 



Ti-so ~ 



a a 
a a 

-52(E<)(E<) (28) 



which is not integrable in general for arbitrary couplings 
(ffi) <?2, 34)- However, in the present case, where ij ^ 1 
ov J ^ U, the effective coupling constants exhibit an 
interesting hierarchy: 



1341 > 1321 > 



(29) 



As a consequence of the above hierarchy, one expects the 
gap in the orbital sector (~ e""'/'^*') to be much larger 
than any other energy scale in the problem. Therefore 
a mean-field decoupling of the interaction term in the 
Hamiltonian H28|) is sensible. We can carry out the mean- 
field decoupling by introducing two Hubbard-Stranovich 
(HS) fields. For that purpose we rewrite the interactions 
terms of Eq. (|28|l as follows: 



^'i(E<) +54(E<) +5<E<)(E<) 

a a a a 

= Ai(Atai -t- A,6i)2 + X2{Ata2 + A.&a)', (30) 



where the notations are 



Ai - \94\ + ^, A2^2|.gi|--f2- 



(ai,6i) = ( 



(02,62) = (- 



154 1 



I52I/2 



I.92I/2 



I54| 



m'+9l x/(f)^+5| 



= )• (31) 



In Eqs. H31(l the hierarchy Eq. (|29|l was employed to sim- 
plify the expressions. Cleary, Ai IA2I. Let us assume 
that A2 is positive. Next two HS fields, ct, C,are intro- 
duced to decouple Eq. ll^ . The resulting Hamiltonian 
in the action form can be written as 



dxdr 



1 



4A2 ^ 



4Ai 

dr - iutdx 
-i{aai + Ca2) 

dr - iusdx 
~i{abi + C62) 



(32) 



i{aai + (a2) 
dr + iutdx 

1(0-61 + C62) 
dr + iusdx 





'pa 




pa 
>tL_ 



pa 
pa 



Now the Majorana fermions can be integrated out ex- 
actly, and the effective action of tr, C is obtained. The 
saddle point approximations! of the effective action of 
CT, C gives 



\92\ /A2.Ut 



(0 - (.>|^(?(^-l) + ^^ln«).(33) 
2|g4| VAi Us TT Us \g2\ ' 

Note that | (ct) | > |(C) | . The saddle point values of A* , Ag 
are 



(At) 
(A.) 



2Ai' 
1.92 1 (a) 



2|<?4| 2Ai 



ut 3 Ai 
1 In 

Us TT Us 



2|54| 



(34) 
(35) 



Note that the factor of 



111 + 3 A^ln M 

Us TT Us 



in Eq. ^ 
< 0. In 



can be either positive or negative since In 
case the factor is positive we have AgAj > 0, while the 
opposite holds for the negative factor. 

Both the saddle point value and the fiuctuation of ( 
[note the factor C^/A2] are very small compared to those 
of a. Thus, due to the hierarchy Eq. (|23l, the HS field 
C can be neglected in the action Eq. (|32|l . The quantum 
fluctuations of a are large since Ai is large. However, 
for the spin sector the fluctuations are suppressed by a 
factor of 61 ~ |g2|/2|(74| owing to the hierarchy [see the 
last line of Eq. (|32|l ]. Thus, we can take the saddle point 
value of a for the spin sector, while the full quantum 
fluctuations should be taken into account for the orbital 
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sector. Namely, a full integration over a is required for 
the orbital sector. Then the effective Hamiltonians in 
both spin and orbital sectors reduce to: 



where Wg = cr6i ~ |(/2|At, and 



a 

(36) 



(37) 

In this limit, the spin excitations consist of a triplet of 
/ree massive Majorana fermions with mass rrig (or equiv- 
alently off-critical Ising models) that span the spin one 
representation of S0(3)spin- This result can be under- 
stood qualitatively as follows. Since J ^ [/, the Hund 
coupling Eq. (j^J dominate in Eq. |^ and the two spin 
1/2 operators, S^i and 8^2 are effectively bound into a 
spin one state. Thus, one expects that in this limit the 
low energy sector of Eq. Q may be identified with that 
of a doped spin-one chain. ^^■^■^ This result will be recov- 
ered in the next section treating the strong coupling limit. 
We also expect that in this limit J ^ i7, all single parti- 
cle excitations will disappear of the spectrum, since the 
strong Hund coupling tend to pair electrons. This is in- 
deed the case : the kinks of the SO (6) GN model, having 
the spin-orbital quantum numbers of the electron, vanish 
and one is left in the lowest energy spectrum with massive 
Majorana fermions, that is to say magnons. The explicit 
calculation of (2 kp component of) spin-spin correlation 
functions using Ising model formalism [ Eqs. (|43I44|| ] 
shows that there is no sharp resonance. This implies 
that the spin excitations are incoherent. 

The orbital sector itself is described by an S0(3)orbitai 
GN model which is integrable. Contrary to what hap- 
pens in the spin sector, there are no stable (Majorana) 
fermions in the excitation spectrum^-'^ . The kink and 
anti-kink states with mass mt ~ e^"*/'^''! exhaust 
the excitation spectrum of SO (3) GN model. 

We see that both A and B regimes differ in their spec- 
tral properties, which is deeply related to the differences 
in the underlying symmetries. Starting in the A phase 
and increasing the value of the Hund coupling J we pre- 
dict that, above a critical value Jc, the gap in the orbital 
sector rrit becomes much larger than the spin gap mg: 
the low energy excitations are exhausted by the spin ex- 
citations. Above Jc fermionic excitations in the orbital 
sector disappear and one is left with solely kinks and 
anti-kinks. This feature is reminiscent of a decoherence 
phenomenon in orbital- like excitations. This prediction 
can be tested numerically. 



2. Charge Sector 

The fact that a gap opens in the spin-orbital sector im- 
mediately suggests the possibility of charge density wave 



(CDW) and superconducting (SC) instabilities. The cor- 
responding order parameters are given by: 



Ocdw = ^ i'^Raa'^Laa + h.C. (38) 



sc,± 



V'i?lTV'L2i ± V'mtV'LU + (^ ^ i), (39) 



where +/— stands for singulet/triplet SC. 

For the discussion of the correlation functions of the 
above order parameters it is most convenient to use the 
Majorana fermion approach.^ To this end let us start 
with the Abelian bosonization formulae: 



(40) 



where 



(41) 



\/27rao 

i'R(L)acr are boson fields satisfying 

[4'B.,aa, 4>L,l3a'] = j/4 5al36a,a 

Kaa- are the Klein factors which enforce Fermi statistics 
A convenient choice of basis for the boson fields is 



R ^ -I- 4>R^ii + 4>f,2] + 't'fai 



,1T 



,1T 



,1T 



-,2T 



-,2T 



-,2T 



(42) 



Refermionization formulae are given by: 







Vs 






•\/27rao 












\/27rao 






risf 


\/2 . 




y/^TTOo 



(43) 



where rja are new Klein factors. The spin and orbital Ma- 
jorana fermion triplets are given by — (C^jC^,^"^) and 
— (C^iC^jC^); respectively. In the Majorana fermion 
basis, both operators Eq. H38|l and Eq. 1)39(1 are non- 
local, while they take local form in terms of the order 
{(Jo) and the disorder (/ia) operators of the six under- 
lying (off-critical) Ising models associated with the six 
Majorana fermions (^^ ,^") — (C^, •■•,C^)- Using the cor- 
respondence: 



^1^2+«0'iCT2, e 



0'2Ml - 
' 0'5M4 - 



1^2(^1, 
- 1^6^3,(44) 
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where </> = (/>l + fpR ^-nd 9 = (p^ — (pji we find: 

Ocdw ^ AiiA^2Ai3M4M5M6 cos Vtt^c 

+ 0"! (72 (T3 (74 0-5 (76 sin Vtt^c- (45) 

Osc,+ ~ e~'^^'^[/ii^2 + i o'ia'2][Ai3A^4M50-6 + (73 0-40-5/16]. 

(46) 

Osc- ~ e~'^^''[mfi2 + i 0-1 0-2] [0-3/14 ^6 - Ai30-4(750-6]- 

(47) 

We are now in a position to discuss the long distance 
properties of the correlation functions of the above order 
parameters. 

Consider first the regime A. The spin-orbital depen- 
dent parts of Eq. (|45I46I47|I . which has scaling dimen- 
sion |, are expressed in terms of products of six order and 
disorder operators aa and ■ These 2^ operators consti- 
tute a basis for the primary operators transforming in the 
spinor representations of SO (6)1 WZW conformal field 
theory. Among the spinorial primary operators, there 
are two S0(6) singlet primary operators : o-iCT2t73 0-4 0-506 
and /J,i/i2/i3/i4/i5/i6. At this point it is worth stressing 
that the enlarged SO (6) symmetry present at low energy 
in the the regime A is different from the original SO (6) 
symmetry of the noninteracting theory. Two symmetries 
are related by the duality transformation Eq. H27|l in the 
orbital sector. Such a transformation interchanges the 
order and the disorder operators: 

Oa ^/ia, a= (4,5,6). (48) 

Consequently only Osc- contains low energy SO (6) sin- 
glets that can take a nonzero average value. Therefore, 
there exists quasi-long range triplet superconducting or- 
der. 

(OL-(^,T)4c,-(0,0))^— ^ , (49) 

In contrast to the triplet superconducting order, both 
CDW and singlet superconductivity have short ranged 
correlations. A similar analysis can be done when J < 
(antiferromagnetic). In this case the model exhibits a 
singlet superconducting instability rather than a triplet 
one. 

In the regime B the situation is different. From Eq. 
(|34I35|I we find that AsAt can be either positive or neg- 
ative. In case of A^At > 0, depending on the sign of At 
we have at the mean field level either: 

(Ma=l,2,3) = (/^a=4,5,6) ^ 0. (50) 

or 

(0-0=1,2,3) = (o-a=4,5,6) 0. (51) 

From Eqs. I|45I46I47|1 . we conclude that a CDW instabil- 
ity is expected in this case. In case of AgAt < we have 
either: 

(^0=1,2,3) = (o-a=4,5,6) 7^ 0. (52) 



or 

(0'a=l,2,3) = (Ma=4,5,6) 7^ 0. (53) 

Then a triplet superconductivity is expected in this case. 

To summarize, we find that at weak coupling the Hund 
perturbation always opens a gap in the spin-orbital sec- 
tor. However, depending on the relative magnitude of U 
and J we may distinguish between two qualitatively dif- 
ferent regimes. At small Hund coupling, J < U, the spin- 
orbital sector displays an effective low energy with en- 
larged SU(4) symmetry and there also exists a triplet su- 
perconducting instability in the charge sector. For large 
Hund coupling, J ^ U, there is a partial SU(2)orbitai 
restoration, and depending on the relative magnitudes of 
parameters the superconducting instability may disap- 
pear in favor of a CDW quasi-long range order. The way 
these two regimes are connected is a nontrivial problem 
and requires a nonperturbative approach. 

IV. STRONG COUPLING 

The situation at large coupling depends on the filling 
as well as the possible umklapp termsj^ For incommensu- 
rate fillings all umklapp operators are strongly oscillating 
and may be discarded at low energy. Consequently, the 
charge excitations are expected to remain massless for all 
U and J. For commensurate filling, i.e. when n = p/q, 
where p and q are co-prime numbers, the possible umk- 
lapp operators allowed by parity and translational invari- 
ance are of the form: 

-ffumklapp = cos[4g\/7i^0c]. (54) 

These operators represent the processes that conserve to- 
tal momentum up to an integer times the Fermi momen- 
tum kp in the continuum limit, and they have the scaling 
dimensions: 

Ag = 4K,q\ (55) 

The equation (|55|l implies that the umklapp operators are 
irrelevant as far as Kc > Xjlq^ . From Eq. H15|l . we find 
this is the case as far as U and J are small enough. How- 
ever, as the coupling constants increase, we expect to 
decrease and possibly to reach the critical value l/(2g^). 
Below the critical value of i^c, the umklapp operator Eq. 

becomes relevant and a gap opens in the charge sec- 
tor. Thus, the mere existence of a Mott-Hubbard (MU) 
transition is related to the nonuniversal dependence of 
K(. on the coupling constants C/, J, and the filling n. The 
dependence is not well-known in general. At present, K^. 
is only known at gwarier-filling and for J = 0.'* In this 
case, Kc reaches its critical value = 1/2 at the value 
V = 2.8t, where a MU transition toward an insulating 
phase has been shown to occur. What happens when 
J 7^ and for other commensurate fillings remains an 
open question. 
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In this section we shall focus on the quarter-Siled case 
and make a reasonable hypothesis that a Mott transition 
still takes place in the presence of a Hund term. Conse- 
quently we expect the model described by ^ to be an 
insulator for large U and J when n — I. For this par- 
ticular filling the strong coupling regime is best achieved 
by going to the Heisenberg limit. As shown by Arovas 
and Auerbach'^'* all of the relevant low energy states at 
strong coupling regimes are given by: 



1 



1 



U2 : c{|c{jO), 4^42 JO) 



(56) 



The above states represent inter-orbital spin-singlet, 
inter-orbital spin-triplet, and intra-orbital spin-singlet 
states. The energy of each state is given by: 

3 1 
uo = U+-J, ui = U- - J, U2 = U. (57) 

The effective strong coupling Hamiltonian depends cru- 
cially on the value of Mi — U — J/2. 



A. [/ > J/2 > t 

In this case all the energies in Eq. H57|l are positive 
and we can employ the results derived by Arovas and 
Auerbachi^ They found that the effective Hamiltonian 
is given by a generalization of the spin-orbital model: 

i 

+ AiS, ■ S,+iT, • T,+i -t- A^TtT^^^S, ■ S.+ij^S) 

where 



and 



Ai 

A2 

A3 

Ai 
As 



2^ (1 - 3??') 
U (l-77)(l + 377)' 
2f_ (1 + 5??) 
IT (l-r;)(l + 377)' 

2t^ 3ri 
^IT (1 + 37;) ' 
2*2 4(1 + r]) 

IT (l-77)(l + 37?)' 

-4^3, 



with 77 = ^ ^ 1. When r] = 0, the Hamiltonian Eq. 
(|58|1 is the SU(4) invariant antiferromagnetic Heisenberg 
model studied by Sutherland^. Expanding the Hamilto- 
nian Eq. (|58|1 with respect to the SU(4) symmetric point 
we obtain 



Ti-cS ~ Ws(7(4) 
2thj 



U 



- 4S, • S,+iT, • T,+i + 12T,'Tl_^,S, ■ S,+iJ . (61) 

Following Ref. j5| we find that in the continuum limit 
the effective Hamiltonian Eq. H61|l is identified with that 
of Eq. H10() . When 77 < 1, a RG analysis reveals the same 
enlargement of SO (6) symmetry as in the A regime of 
the weak coupling case. Thus, one may conclude that the 
SO (6) symmetric A regime extends from weak to strong 
couplings for small enough Hund interaction. 



B. 



t < (7 < J/2 



In this case ui becomes negative and the strong cou- 
pling approach developed by Arovas and Auerbach does 
not apply. However, another strong coupling expansion is 
sensible when [/ <C J/2. Indeed, in this limit the ground 
state consists of the local spin triplets which contain two 
electrons per site and spontaneously break translational 
invariance. Assuming that the spin triplets are located on 
even sites (such that odd sites are empty) we can find an 
effective interaction between the local spin triplets. The 
effective interaction can be determined from the strong 
coupling expansion to the order of i"* in a straightfor- 
ward way. Let us denote the spin triplet at the 2j-th 
site by I2J, which is a S" = 1 spin operator. By solving 
the associated two-site problem up to the order of t"*, the 
following effective Hamiltonian can be obtained: 



'T-cff — 



(62) 





a 

.2". 


(T(T' 


J 

with 


/ y ^laa- 
craa.' 


"cr" 
.2"- 


aa' 


(59) K, = 



7-31e/2-K2e2 + e'V2 



5J3 (1 
24*4 



)3(l-e/5)(l-eV4)' 
1 



5J3 (l-e/5)(l-e)2 



(63) 



(60) 



where e = ^ < 1. The Hamiltonian Eq. is 
that of the antiferromagnetic Heisenberg spin 1 model 
with exchange Ki and a biquadratic coupling K2 < 0. 
In this regime of coupling constants, the excitations are 
massive and consist of a triplet of Majorana fermion^S^ 
It is however hazardous to conclude that the Haldane 
magnons constitute the low energy excitations in this 
strong Hund coupling regime. Indeed, the groundstate 
for the t = 0,U <^ J limit of the model is doubly degen- 
erate for the local spin triplets can be located on even or 
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odd sites. In addition to the Haldane magnons of H62() 
we thus have to take into account the kinks that connect 
these two groundstates. Note that these kinks carry in- 
teger spin, since excitations built out of single electrons 
have very high energy in this C/ <C J regime. 

A rough estimate of the energy of kink excitations can 
be obtained in a static approximation. Consider a pined 
kink-antikink pair on top of the groundstate of H62f) . lo- 
cated at sites p and q {p, q beeing odd integers) , that is 
to say that the local spin triplets are located on even 
sites 2i for 2i < p and 2i > q and on odd sites 2i + 1 
ioi p < 2i + 1 < q. When the t perturbation is in- 
cluded, this will result in the effective hamiltonian H62() 
with modified exchange Ki^2 on the bounds p and q. The 
static approximation is thus equivalent to two bound de- 
fects, one weak and one strong (with transparent nota- 
tions, 1^1.2 (p) ^ i^i,2 ^ -^1,2 (<z) for the particular con- 
figuration we have chosen). The effect of bound defects 
has been studied in the spin 1 pure Heisenbcrg chain 
(without biquadratic exchange), which has the same low 
energy physics as (1^ . The defects there lead to the ap- 
parition of a triplet excitation inside the Haldane gap^ 
It is more than likely that beyond this static approxi- 
mation, the kinks will form a band of spin 1 excitations 
once kinetic energy is included, as occurs for holes in the 
AKLT model.27 

While this problem would clearly require a more care- 
ful analysis, which goes well beyond the scope of this 
paper, this simple static picture indicates that the lowest 
energy excitations will be massive spin 1 kinks. In any 
case, the low energy spectrum is exhausted by massive 
magnons. The effective continuum theory in regime B at 
weak coupling t ^ J ^ U describes also massive spin 1 
particles ; it is thus tempting to postulate a continuity 
from weak coupling to strong coupling for the low energy 
excitations. 



C. t/ ~ J/2 » i 

In this case the local spin triplet two particle states 
have very low energy and they will mix with other states 
with one particle per site in the first order of hopping t. 
These states should be diagonalized first before taking 
into account the higher order perturbation in t. The 
detailed study of this regime is beyond the scope of this 
paper. 

V. SUMMARY 

We have studied the one-dimensional SU(4) Hubbard 
model away from half-filling perturbed by a Hund cou- 
pling J in both the weak and strong coupling limits. The 
Hund coupling turns out to be relevant irrespective of the 
short ranged Coulomb repulsion U and is responsible for 



the opening of a spectral gap in the spin-orbital sector. 
We found two qualitatively different regimes depending 
on the relative strength between J and U. 

For small enough Hund coupling J, we found that the 
symmetry is dynamically enlarged to SU(4) at low en- 
ergy. In this limit, the the spin-orbital degrees of freedom 
are described the integrable SO (6) Gross-Neveu model 
and the excitations are found to be coherent. In this 
regime, we found a superconducting instability when the 
charge excitations are not gapped. For large Hund cou- 
pling the SU(4) symmetry is no longer enlarged. Instead, 
we find that the orbital degrees of freedom decouples at 
high energy and the low energy excitations lie in the spin 
sector and consist of the three free massive Majorana 
fermions. As a result, the spin excitations become inco- 
herent. Our results are summarized in Fig. ^ 
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FIG. 1: The phase diagram for the spin-orbital degrees of 
freedom. The symmetry as well as the effective model at the 
strong coupling (low energy) regime are indicated. Bound- 
aries represent smooth crossover rather than critical quantum 
phase transitions. 
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